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Abstract 

Dark matter, X, may be generated by new physics at the TeV scale during an early matter- 
dominated (MD) era that ends at temperature Tr <C TeV. Compared to the conventional 
radiation-dominated (RD) results, yields from both Freeze-Out and Freeze-In processes are 
greatly suppressed by dilution from entropy production, making Freeze-Out less plausible 
while allowing successful Freeze-In with a much larger coupling strength. Freeze-In is 
typically dominated by the decay of a particle B of the thermal bath, B —> X. For a large 
fraction of the relevant cosmological parameter space, the decay rate required to produce 
the observed dark matter abundance leads to displaced signals at LHC and future colliders, 
for any mx in the range keV < mx < m R and for values of ms accessible to these colliders. 
This result applies whether the early MD era arises after conventional inflation, when Tr 
is the usual reheat temperature, or is a generic MD era with an alternative origin. In 
the former case, if mx is sufficiently large to be measured from kinematics, the reheat 
temperature Tr can be extracted. Our result is independent of the particular particle 
physics implementation of B —>• X, and can occur via any operator of dimension less than 
8 (4) for a post-inflation (general MD) cosmology. An interesting example is provided by 
DFS axion theories with TeV-scale supersymmetry and axino dark matter of mass GeV to 
TeV, which is typically overproduced in a conventional RD cosmology. If B is the higgsino, 
h, Higgs, W and Z particles appear at the displaced decays, h —»• ha,Za and /A —>• W^a. 
The scale of axion physics, /, is predicted to be in the range (3 x 10 s — 10 12 ) GeV and, 
over much of this range, can be extracted from the decay length. 
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1 Introduction 

In the early universe, when the background radiation had a temperature from eV to MeV, 
the energy density was dominated by radiation. How far back did this radiation dominated 
(RD) era extend? When inflation ended, provided the inflaton decays were not extremely rapid, 
there was an era of matter domination (MD) which ended at the reheat temperature Tr. It is 
commonly assumed that Tr is very high, many orders of magnitude above the TeV scale, but 
observationally the most stringent constraints are from the effective number of neutrino species 
from CMB measurements, Tr > 4 MeV [1], and from Big Bang Nucleosynthesis [2]. This early 
MD era could extend for many decades in temperature above T R , and include the TeV epoch. 
There are several alternative origins for a long early MD era, including long-lived heavy particles 
that were once in thermal equilibrium and oscillating fields composed of light bosons. When 
the MD era results from inflation it has an evolution that is purely Non-Adiabatic in character, 
giving a MD^ era, as shown in Fig. 3. The more general MD era splits into two, beginning 
with Adiabatic evolution, MD^, and ending with MD^a, as illustrated in Fig. 1. 

In this paper we explore the consequences of having the dark matter abundance determined 
by reactions during an early MD era. This profoundly changes the abundance compared to the 
conventional assumption of dark matter genesis during the RD era for two reasons. Firstly the 
time-temperature relation is changed so that any given temperature occurs at earlier times and 
secondly, once the genesis reactions stop, the comoving abundance of dark matter is diluted by 
entropy production. The second effect dominates, so that the MD era results in a depletion 
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of dark matter relative to the usual RD result; by a factor (T R /T F ) 3 for Freeze-Out during 
MDjv A [ 3, 4, 5], where Tp is the Freeze-Out temperature, and by a factor T R /(TmTp) 1 ^ 2 for 
Freeze-Out during MD^ [6] where Tm is the temperature at which MDa begins. This factor is 
likely very small, suggesting a dark matter mass much larger than the TeV scale. In this paper 
we focus on TeV scale physics as the origin of dark matter, as this could be probed at colliders, 
and therefore we seek an alternative to Freeze-Out, as Freeze-Out during a MD era is unlikely 
to yield sufficient dark matter. 

Particle dark matter requires an addition to the Standard Model: a cosmologically stable 
particle X of mass nix ■ How is the abundance of X determined? The Freeze-Out mechanism 
results if X has sufficient interactions with the known particles that it is in thermal equilibrium 
at temperatures T of order rrix■ As T drops below mx, X tracks a Boltzmann distribution for 
a while but, as it becomes more dilute, its annihilation rate drops below the expansion rate and 
it freezes out of thermal equilibrium. This mechanism has great generality, applying to a wide 
range of theories where the interactions of X are sufficient to put it in thermal equilibrium at 
T ~ mx- Furthermore, it is highly predictive since it is an IR dominated process - it does not 
depend on any physics at energy scales above its mass. Indeed, Freeze-Out during the RD era 
suggests mx is broadly of order a TeV and that X has an interaction rate that can lead to 
signals at both collider detectors and direct and indirect detection experiments. 

Freeze-In provides another general production mechanism and results when the interactions 
of X are insufficient to bring it into thermal equilibrium [7]. 1 At temperatures T above mx 
these feeble interactions allow the production of X from decays or scatterings of some bath 
particle B of mass m R at rate r(T). This produces a yield of X particles at time t{T ) of 

Y% rod (T) ~ r(T) t(T) (1) 

which is IR dominated when the interactions between B and X are of dimension 4 or less. 
Taking B to be the lightest observable sector particle carrying the stabilizing symmetry, the 
production rate becomes Boltzmann suppressed below mg, so that the dominant contribution 
to Y x arises at T ~mg. For T > mg, the yield Y(X) grows towards its equilibrium value, but 
it never reaches equilibrium, and this Freeze-In towards equilibrium stops once T drops below 
mp- Decays generally dominate over scatterings, so that in this paper we study the reaction 

B —>■ AsmA" (2) 

where Ag M is one or more Standard Model particles. 

When Freeze-In occurs during the RD era, the value of the B lifetime that leads to the 
observed dark matter abundance is of order 10 _2 s for mp,x of order the weak scale. This leads 

to interesting signals of long-lived stopped particles at colliders if B carries electric or color 

charge [7], but not if B is neutral. The interaction between B and X is described by a very 
small coupling of order 1CT 12 , as shown by the blue curve of Fig. 4. 

In this paper we study Freeze-In during an early MD era, requiring T R to be less than mg, 
and fold that the B lifetime, tb, is reduced. Freeze-In via scattering during the MDa^ era was 

: We do not consider the case of dark matter arising from decay of the matter that dominates the energy 
density during the MD era [8]. 
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studied in [4, 5], whereas Freeze-in via decays for goldstino DM was studied in Ref. [9]. The 
reduction in tb arises because the relevant time scale t(ms) that appears in (1) is smaller in 
the MD era than in the RD era and because the produced yield Y£ rod (T ~ms) is subsequently 
diluted by entropy production. Both effects require an increase of the rate T appearing in (1). 
In a model-independent analysis, we find that the required tb depends on two particle physics 
parameters () and two cosmological parameters (_D,Tr), where D appears only if the 
MD era has a MD^ portion, in which case it represents dilution. For Freeze-In during a MD^ 
era we show our results in the two panels of Fig. 6 for values of m b relevant for LHC and 
future hadron colliders, respectively. The decay length tb in this case is function of the reheat 
temperature T R only. If we also allow a MD^ era then we have contours for Tb in the (D, T r ) 
plane, as show in Figs. 8 and 9. Over a large portion of the parameter space for Freeze-In during 
the MD era via (2), these B decays can be probed at colliders as they lead to exotic signals 
displaced from the primary vertex. This result is remarkably general, applying to a wide range 
of dark matter masses, keV < mx < ms, and a wide range of interactions that induce the decay 
(2) 2 . For Freeze-In during the MD^ era, any interaction of dimension 4 or less is sufficient, while 
the extreme UV-insensitivity of the MDjva era allows any interaction of dimension less than 8. 3 
While the interaction is still feeble, the strength is several orders of magnitude larger than in 
the RD case; for example, a coupling of order 10“' for the dimension 4 case, as in the red curve 
of Fig. 5. 

The displaced collider signals have a size depending on the production rate of B and a 
character depending on the identity of A S m and tb- For B color and charge neutral, possibilities 
for Asm include (h, Z, d + i~, qq, 7 ) giving a variety of non-prompt decays. Using Run I data, LHC 
experiments have placed limits on displaced vertices [11, 12], displaced jets [13, 14] and displaced 
leptons [15, 16]. However, some of these analyses have assumed particular forms for the signal, 
and may not be sensitive to, or optimized for, Freeze-In. For example, it would be interesting 
to search for displaced signals involving a Higgs or Z boson, as arise from B —y hX, ZX. For 
B carrying color and/or charge, Ag M includes {W, £ ± , q, tv, qq, g) giving decays from long-lived 
particles stopped in the detector [17, 18]. 

Axino dark matter in supersymmetric DFS theories of the QCD axion provides an interesting 
example of Freeze-In. With Freeze-In during the conventional RD era too much dark matter 
is generically produced, unless theories are crafted to have very high axion scales, /, or very 
low axino masses. However, for Freeze-In during the MD era, axinos account for the observed 
dark matter for / in the range (3 x 10 8 — 10 12 ) GeV and axino masses in the GeV to TeV 
range. Furthermore, as shown in Figs. (10) and (11), over a large fraction of the parameter 
space not excluded by axion overproduction and white dwarf constraints, there are displaced 
signals at colliders that arises from the higgsino decays that induce the cosmological abundance, 
h —» ha, Za and hA —» W ± a. 

2 Our results also apply to B —> XA', with A' containing exotic non Standard Model particles, followed 
by rapid decays A' —> A SM . The Freeze-In yield and the collider decay length are then both determined by 

r(R -> xx). 

3 Our results also apply if X is part of a hidden sector and FI dominates dark matter production [10]. 


4 



2 An Early Matter Dominated Era 


In this section we review the evolution of the universe during an early Matter Dominated era, 
MD, and study the parameter space of such a cosmology relevant for Freeze-In (FI) of dark 
matter. We consider an FRW cosmology with two components: Radiation (R) and unstable 
Matter (M), that decays to radiation with a decay rate Tm- The radiation is thermal at 
temperature T and has energy density 

(3) 

where g = g{T) counts the effective number of relativistic species at T, including particles 
from the Standard Model and beyond. The matter M can be fermonic or bosonic and, in the 
era of interest to us, is not thermally coupled to the radiation. The M particles might have 
been in thermal equilibrium at an earlier epoch when they were relativistic, so that they have a 
number density comparable to that of the relativistic species, or they might arise from a classical 
oscillating bosonic field, as with inflatons or moduli, in which case their number densities could 
be much larger. When T is of order the MeV scale M have decayed and the universe is in the 
conventional Radiation Dominated era, RD, but at earlier times we insist that the initial density 
of M is sufficiently large that there is a long matter dominated era. 

The Boltzmann equations governing the evolution of the energy densities Pr,m are 


dpM 

dt 


+ 3 H pm 


dpn 

dt 


+ 4 Hp R 


“ F M pM 


T m Pm 


( 4 ) 


where H is the Hubble parameter. 4 The solutions to these Boltzmann equations depend on 
the initial densities Pr^m, and Tm- At times earlier than tm = 1/Tm, the effect of decays on 
Pm can be ignored, giving 


Pm 

Pr 





1/2 

F mM P i p 1 ^ 




( 5 ) 


where a is the Robertson-Walker scale factor and Mpi = G -1 / 2 = 1.22 x 10 19 GeV. There are 
two contributions to the radiation energy density, the first is the red-shifted initial density while 
the second, proportional to T M , arises from M decay. 

This solution is illustrated in Figure 1 for a* sufficiently small that the universe starts in an 
early radiation dominated era, RD'. It transitions at to the early MD era, which ends at 

4 The second line of Eq. (4) is derived from the first law of thermodynamics under the assumption that the 
equation of state pp = pr/2> is satisfied. Strictly speaking, this is not exactly true when the temperature is near 
the mass of any of the particles in the thermal bath. We verified that the corrections to the equation of state 
are small and they only affect the final results by a few percent. Therefore, we take a constant value of g * in our 
numerical study. Freeze-in is dominated at temperatures of the order of the weak scale, so we take g = 106.75, 
accounting for the full Standard Model degrees of freedom. 
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Figure 1: Cosmological history with a, < a^. A long early MD era commences at clm] it 
has an initial regime with adiabatic evolution and then at a^A it transitions to Non-Adiabatic 
evolution, where Pr is dominated by decays of M even though t <C Tm • Finally, the majority of 
M decay at an where the conventional RD era begins. 


reheat, or, when M decay and pM drops to zero. The MD era is divided into two regimes. In the 
earliest one pr is dominated by red-shifted initial radiation, so that the second term in pr can 
be ignored and the evolution is Adiabatic (the MD^ era). The second regime has pr dominated 
by radiation from M decays, and hence the evolution is Non-Adiabatic (the MD^ era); the first 
term in p R can be ignored so pr ~ TmMpi Pm ■ F PRi i s sufficiently large the universe evolves 
through both MD regimes; otherwise there is only a MD^ regime. The transition between 
these two types of cosmology occurs at 

PRi ~ (TmMri) p)£. ( 6 ) 


2.1 General MD Histories Relevant for Freeze-In 

In this paper we study Freeze-In of dark matter during this MD era. Taking ^ matches 

on smoothly to the RD' era, but it could be that the Boltzmann equations of (5) break down 
at some a larger than cim, so that the MD era emerges from a non-RD era, such as inflation or 
topological defect domination. To be completely general we allow a* to be anywhere between 
(im and cir. (Taking a* -C aw gives an identical MD era to the case of a,i = om, and cq > cir 
does not have an early MD era.) Equivalently, we study all initial conditions with 

PMi > PRi, (FmMpi) 2 . (7) 

In fact we insist that pu % S> (T m^pi ) 2 so that the MD era is not a transient phenomenon. 

The set of cosmological histories of interest to us is shown shaded in Figure 2, and splits 
into two according to whether the evolution is initially MD^ (light gray) or MDjva (dark gray). 
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Figure 2: The space of cosmological histories in the (prt n Pm { ) plane. The unshaded region below 
the dashed line has no early MD era. The three regions above the dashed line are distinguished 
by the initial era. The blue and red trajectories give the evolution of energy densities for 
histories labelled A and B (with axes interpreted as Pr(cl) and Pm(u), a > <R, and several key 
temperatures indicated). Note that T^Mj,; ~ T^. 

The region above the dashed line which is initially RD' will evolve to give the same MD era as 
histories on the boundary between initially RD' and initially MD^. Two particular histories are 
illustrated by the red and blue trajectories having initial energy densities A and B. Point A has 
Pr i = pMi (r a jMpi) 2 and corresponds to cq = czm in Figure 1, leading to both Adiabatic and 
Non-Adiabatic evolution. Point B has pR t = 0. In this case, those M that decay much earlier 
than tm rapidly increase the radiation density to pr ~ T m M P i p^ , leading to a MD era that is 
purely Non-Adiabatic. The horizontal part of the blue trajectory is covered very quickly so that 
the MDTvn evolution occurs on the non-horizontal part of the blue trajectory. For fixed Tm, the 
possible histories with purely Non-Adiabatic evolution can be described by pMi, which labels 
points along the line of non-adiabatic evolution. Histories having both evolutions are described 
by ( PMiiPRi ), corresponding to points in the region with light shading. 

Histories with pr. = 0, such as the blue trajectory in Fig. 2, could emerge from an inflationary 
era with M interpreted as the inflaton, which starts to oscillate at cq, as illustrated in Fig. 3. 

2.2 Evolution of time, temperature and energy density 

For the purposes of studying Freeze-In during the MD era, we lose no generality in considering 
histories with p M . = p R . (corresponding to a* = a^, such as point A in Figure 2) since, in the 
late part of their evolution, they include all other histories that have larger a*. In particular, 
solutions with equal initial conditions always contain a MDjvn phase at late times; therefore they 
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Figure 3: A long MD era after inflation. 


also cover histories with only a MDjva era (up to the very quick beginning of it, shown on the 
horizontal blue line in Fig. (2)). Hence the three dimensional parameter space (pn v pu. t ,Tm) 
may be collapsed to two (pa., = The evolution of time and temperature in these 

histories, relevant for physical processes, is discussed below and the two parameters may be 
taken as (Tm, Tr), the temperatures at the beginning and end of the MD era. They are related 
to (pr^Tm) by p ~ TfIf and Tr ~ V^mMpp 

During the entire MD era the age of the universe is t — tM{ci/ clm ) 3 ^ 2 ■ During the MD^ 
regime g l ^aT is constant. For a simple pedagogical discussion, we will assume a constant g so 
that the energy density is 

pa(t) = (r 4 + r„r 3 ) (8) 

and in this regime the time-temperature relation has the parametric form 


t A (T) ~ 


1 Mpi 

g 1 / 2 T ^ 2 T 3 / 2 


( 9 ) 


The universe transitions from adiabatic to non-adiabatic evolution when the two terms in 
Pr of (5) are comparable, at time 5 

t ,2/5 ,3/5 nn x 

tNA ~ t R (1UJ 

with t M = t A (T M ) and tR = tm = T^ 1 , and at temperature 


Tna ~ T m 



( 11 ) 


5 The ~ symbol indicates approximations resulting from complicated behaviors during transition regions be¬ 
tween eras, while the ~ symbol indicates that, in addition, numerical constants have been dropped. 














In Fig. 2, points along the boundary separating the initially MD^ and MD^a regions can be 
labeled by Tma- For histories resulting from inflation there is no MD^ era, rather the history 
begins at T^a, which is identified with 

Tmax ~ VEjTji ( 12 ) 

where Ej is the energy scale of inflation. For this case it is convenient to take the independent 
variable as T max instead of Tm- 

During the MDjva regime a 3 / 8 T is constant so that the energy density is 

= + (13) 

where p\j has been normalized so that it is equal to p R at T R . This gives a time temperature 
relation 

, m MpiTr n4 x 

tNA\T) ~ gl/2rp4,' 

Matching (8) and (13) gives 

t%a ^ r„T*. (15) 

The final RD era is reached when T = T R and t = t R . Including a numerical factor, (14) at 
this transition gives 



Above T eq ~ leV the universe was radiation dominated and it is frequently assumed that 
this RD era began at extremely high temperatures, as high as the energy scale of conventional 
inflation, Ej, which could be as large as 10 16 GeV. However, in the cosmology discussed here, 
the universe had a long MD era from T R to Tm- Observational constraints are very mild: T R > 3 
MeV from CMB and BBN and Tm Ej. Even in the case that this early MD era arose from 
the conventional inflaton, and hence was purely MDtva, it could have lasted for many decades in 
temperature from T R as low as 3 MeV to T max ~ \JEjT r . If the DM abundance was determined 
during the MD era we expect significant changes in predictions compared to the RD case for 
two reasons. The DM yield depends on the Hubble parameter at production (Y oc H for FO 
and Y oc 1/H for FI) and at any temperature in the MD era H is larger than at the same 
temperature in the RD era by an amount that can be as large as (Tm/Tr) 2 ^ when T = T^a- 
However this is overwhelmed by dilution effects subsequent to DM production, which reduces 
the yield by as much as Tr/Tm■ For the production of TeV-scale dark matter in a long MD era 
this is typically catastrophic for FO since the abundance becomes too low, but for FI it has the 
interesting feature of enhancing the required couplings of the weakly interacting state to yield 
displaced signals at colliders. 
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3 Freeze-In During an Early MD Era 


We consider Freeze-In (FI) production of dark matter during the early MD era via the decay 
B —> AsmW. B is the lightest observable sector particle that carries the stabilizing symmetry, X 
is the dark matter and Asm is one or more Standard Model particles. In addition to FI, there is 
a Freeze-Out (FO) population of B that eventually decays to X. In this Section we focus on FI, 
deferring FO to Appendix B where we show that in the early MD era it is typically negligible. 

The number density evolution of X is described by the Boltzmann equation [7] 


dn x 

dtt 


+ 3Hn x = T B n e £ 


Ki[m B /T] 
K 2 [m B /T} ’ 


(17) 


with Tb the width of B and Ki j2 [x] the first and second modified Bessel functions of the 2nd 
kind. The equilibrium number density obtained using Maxwell-Boltzmann statistics reads 


n 


eq _ 
B — 


9b 

2tt 2 


m 2 B T K 2 [m B IT] . 


(18) 


At high temperatures (T m B ) we recover the T 3 abundance for a relativistic species, whereas 
at low temperatures (T C mg) the number density has the Maxwell-Boltzmann exponential 
suppression. For this reason, the FI production of X is dominated at temperatures Tpi ~ m B . 

If FI occurs during the RD era, the Boltzmann equation (17) can be easily solved [7] giving 
a final yield for A" 


n x 


Y x = — = 4.4 x 10~ 12 ( — 


g B \ ( 106.75\ 3/2 /300 GeV\ ( V B /m B 


9 * 


m B 




1.8 x 10~ 25 


(19) 


where s is the entropy density. Observations fix the DM density but the yield is fixed once the 
mass of X is known using 


£,dm — -— ixidmYdm — 0.44eV (20) 

s 

which is close to the usual temperature of matter radiation equality, T eq ~ leV. We can thus 
rewrite eq. (19) as 


t _ c (9 b\ ( 106.75\ 3 ^ 2 / rri x \ /300 GeV\ f T B /m B \ 

Cy - Zdm [ 2 ) ^ gjf J v!00 GeV/ { m B J V 1 - 8 x 10 " 25 / ' 


( 21 ) 


The coupling A, defined by 

A 2 

Tb = ~^~ m Bi (22) 

8tt 

must be very small to avoid overclosure. For the reference masses and number of spin states 
shown in (21) the observed DM density results for A ~ 2 x 10 -12 . 

In Fig. (4) we show the FI solutions for the DM yield Yx(^)> with time variable x = m B /T. 
We show three different lines corresponding to three different values of A (namely three different 
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m B = 300 GeV m x = 100 GeV g B = 2 



Figure 4: Yields Yx(x) for FI during the RD era for three values of the decay width, as defined 
in Eq. (22), with fixed reference values m B = 300 GeV, nix = 100 GeV, and g B = 2. 

decay widths as in Eq. (22)). In all cases we see that the process is maximally active at 
temperatures of order m B . For the reference parameters of (19), the blue line reproduces the 
observed DM density, whereas the brown and red lines overproduce the amount of dark matter. 
In particular, the red line reaches an asymptotic value of Y\ which is 10% the equilibrium value 
V£ q , and above this point the approximation that DM arises from FI through the decay in 
Eq. (2), and without the inverse reaction, breaks down. 

The FI abundance depends on the cosmological evolution at and after FI. The solution of 
the Boltzmann equation (17) is given in Appendix A for an arbitrary cosmological evolution in 
Eqs. (58) and (59). In this section we apply this to the cosmological backgrounds discussed in 
the previous section, with energy density p R + pm evolving according to the Boltzmann equation 
(4). These cosmological backgrounds depend in general on three parameters, prt t , pM t and Tm 
(or equivalently T R ). In the next subsection we study the case of FI during the MD^ era, 
which is particularly simple because there is no dependence on p R . or p M . and the results can 
be displayed as a function of T R . 

3.1 Freeze-In during Non-Adiabatic evolution 

Freeze-In during MD^ occurs between a^A and a R , as shown in Fig. 1, and the final DM 
density depends on only one parameter of the background, the reheat temperature T R . This 
includes the important example of FI during reheating after inflation, as sketched in Fig. 3. In 
this case inflation could be the conventional one with more than 60 e-foldings with the matter 
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m B = 300 GeV m x = 100 GeV g B = 2 



Figure 5: Yield Yxix) for Freeze-In during the MDjvyi era, for example during reheating after 
inflation. For each of three A, Tr is chosen to give the observed DM abundance, as shown in 
Eq. (23). We take fixed reference values mg = 300 GeV, mx = 100 GeV, and gB = 2. 

M identified as the usual inflaton, or it could be a later era of inflation with fewer e-foldings. 

The FI yields Yx(x) are shown in Fig. 5 for three different values of the decay width 

f 2 x 10“ 12 T r = 1 TeV 

A = < 2 x 10- 8 T r = 4.6 GeV . (23) 

{ 2 x 10“ 7 Tr = 2.4 GeV 

In each case we choose Tr in such a way that we reproduce the observed DM abundance. The 

values of A span a wider range than the case of Fig. 4 since Freeze-In during a MD era allows 

for larger couplings. The blue line has Tr > iur, so that FI occurs during the RD era. As A is 
increased, the FI process becomes more powerful, as illustrated by the brown and red lines, and 
to obtain the observed DM abundance much lower values of Tr must been taken, so that the A" 
abundance is diluted by entropy production after FI. However, it is not possible to arbitrarily 
increase A and still be in the Freeze-In regime. For sufficiently large A, the peak of the Y x 
functions in Fig. 5 will reach the equilibrium value and FI is no longer the applicable production 
mechanism. For reference values {rn B — 300 GeV, mx = 100 GeV, gB = 2) this gives a lower 
limit Tr > 0.6 GeV. At the same time we have to keep Tr < mg; otherwise FI occurs during 
the RD era and there would be no dilution. 

The B decay length necessary to reproduce the observed DM abundance is shown in Fig. 6 
as a function of Tr. The left (right) panel is for a value of = 300 GeV (3 TeV) relevant for 
the LHC (future colliders) and four curves are shown in each, giving results for a wide range 


12 







m B = 300 GeV g B = 2 



0.0 0.5 1.0 1.5 2.0 2.5 

l°9io r fi( GeV ) 


m B = 3 TeV g B = 2 



1.0 1.5 2.0 2.5 3.0 3.5 4.0 

l°g 10 7fi(GeV) 


Figure 6: B decay length as a function of Tr for FI during the MD^ta era, as occurs after 
inflation, for two different choices of m R - In each panel we consider four different values of mx- 
The various curves stop on the left because Freeze-In no longer occurs as A" thermalizes and on 
the right because Freeze-In occurs during the regular RD era. Displaced collider signals occur 
in the shaded horizontal bands as described in Table 1. 


of mx- Signal regions for various displaced events at colliders are given by shaded horizontal 
bands, as described in Table 1. The stopped particle decay signal is absent if B is color and 
electric charge neutral. 


Shaded region 

Decay length 

Signature from LOSP 

Neutral 

Charged 

Dark green 

10 2 cm < tr < 10 2 cm 

Displaced vertices 

/ 

/ 

Light green 

10 2 cm < t b < 10 4 cm 

Displaced jets/leptons 

/ 

/ 

Light blue 

10 4 cm < tr 

Stopped particle decays 

X 

/ 


Table 1: Displaced Collider Signals 


3.2 Freeze-In during Adiabatic or Non-Adiabatic evolution 

We extend the previous discussion to the more general framework of Fig. 1, which also allows for 
a period of adiabatic evolution during the MD era. The full parameter space of the cosmological 
background is (p R .,p M ., T M ) but, as discussed in Section 2, for FI there is no loss of generality 
in taking p R . = p Mi . The parameter space may then be taken as ( T M ,T R ), where T M is the 
temperature at the beginning of the MD era. We trade Tm for another parameter, the dilution 
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m B = 300 GeV m x = 100GeV fife = 2 



Figure 7: DM yield Y x = n x /s as a function of x = m B /T for Freeze-In during an early matter 
dominated era. For each of the three A, we choose (Tr, D) to give the observed dark matter as 
shown in Eq. (24). 


D = Sf/Si ~ Tu/Tri where S % and Sf are the entropy per comoving volume before and after 
the matter decay. 

We consider the same three values of the coupling A as in Eq. (23), and we choose the dilution 
factor and the reheat temperature to reproduce the observed DM abundance. There is some 
degeneracy in the process, and for the purpose of illustration we fix 

f 2 x 10" 12 T r = 1 TeV D = 1 

A = < 2 x 10" 8 T r = 4.6 GeV D = 10 8 . (24) 

[ 2 x 1(T 7 T r = 100 MeV D = 4.7 x 10 7 

The solutions are shown in Fig. 7. This time we notice a different behavior between the brown 
and red lines. The reason is a consequence of our choice of (Tr, D ): for the (brown) red line the 
Freeze-In happens during the (non-)adiabatic matter dominated era. 

Figs. 8 and 9 provide a complete exploration of FI in the (Tr, D) plane, with results shown 
for different (m B ,m x ). The LOSP mass values of m B = 300 GeV (Fig. 8) and 3 TeV (Fig. 9) 
are chosen to be relevant for LHC and a future collider, while values of m x are chosen to span 
the very wide range of possible values that is consistent with FI. Remarkably, the displaced 
vertex region, shaded in green, is always within the FI domain and covers a substantial part of 
the parameter space where the observed DM results from FI. 
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Log 10 r e (cm): m B = 300 GeV m x = 100GeV g B = 2 


Log 10 r e (cm): m B = 300 GeV m x = 1 GeV g B = 2 



logio^R(GeV) 

Log 10 r s (cm): m B = 300 GeV m x = 100 MeV g 6 = 2 




Log 10 r s (cm): m B = 300 GeV m x = 1 MeV g e = 2 



logio^RfCeV) 


Figure 8: Contours of the B lifetime tb (in cm) that give the observed DM abundance via FI 
in the (' Tr,D ) plane. We fix mg = 300 GeV and in each panel we consider different values of 
mx- The upper-left gray region does not give FI, as X thermalizes. The three dashed red lines 
separate four regions where FI occurs during (left to right) RD', MD^, MDtva, and RD eras. 
The red and brown stars in the top-left panel identify the benchmark points chosen in Eq. (24) 
and Fig. (7). Displaced collider signals occur in the green and blue shaded bands, as described 
in Table 1. 
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Log 10 r s (cm): m B = 3TeV m x = 1 TeV g B = 2 



I°9io7r(G@V) 


Log 10 r s (cm): m e = 3TeV /r?x = 10GeV g B = 2 



iog 10 T R (GeV) 


Log 10 r s (cm): m B = 3 TeV m x = 1 GeV g B = 2 



logioT'RfGeV) 


Log 10 r s (cm): m s = 3 TeV m x = 10MeV g B = 2 



logi 0 T R (GeV) 


Figure 9: Same as Fig. (8) but for mg = 3 TeV. 
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3.3 Analytic Approximation to Freeze-In Abundance 

Many features of the numerical results presented in the previous two sub-sections can be under¬ 
stood from simple analytic approximations. The FI abundance of X generated during a Hubble 
doubling temperature at T during era i is 

Y Prod (T) ~ r(T) ti(T) (25) 


where Y (T) is the total rate of X production at temperature T. For clarity numerical factors 
are dropped in this sub-section, unless they are very large. The age of the universe at T is 


U(T) 


H,(T) 



x 1 / 2 


T, 


1/2 • 


M 


(26) 


with i — 1—4 running over the (RD, MDjva, MD (, RD') eras, with RD the usual RD era at 
T < Tj; and RD' the early one at T > Tm ■ When FI is induced by an operator of dimension 4, 
the rate of X production is given by 


r (T) 



(27) 


at T > ms- The upper entry of the vector is for decays, with thr/T the Lorentz boost factor, 
while the lower entry is for scatterings. 

This analytic yield Y i Prod (T) approximates the solution to the Boltmann equations, Yx(T ), 
for T > rriB, where dilution does not have a large effect. For the case of decays, the power laws 
in the various eras are Y Prod (T ) oc (1/T 3 ,1/T 5 ,1/T 5 / 2 ,1/T 3 ). At x <C 1, the RD, MD^ and 
MD a slopes are visible in Figs 4, 5 and 7, respectively. (In Fig. 5 the blue curve starts with 
MDat A behavior, but the slope decreases as x approaches unity where FI occurs during the RD 
era. In Fig. 7 the blue curve has D = 1, so that there is no MD era and the slope reflects RD.) 

The contribution to the DM yield today from FI production when the universe had temper¬ 
ature T and was in era i is 

Y'Prod (X'j 


Y„,(T) 


ACT) 


(28) 


DAT) 


(29) 


where Di(T) is the dilution factor from non-adiabatic decays of M during the MDat A era between 
T and Tr 

/ T 1 5 'T' rp 

i z_ 

’ Tf Tr' Tr 

During the MDtva era, Eq. (13) shows that rix oc T 8 , so that the dilution factor of Yx scales 
as T 5 . As the temperature increases to Txa the dilution reaches T^ A /Tff which by Eq. (15) is 
Tm/Tr, the constant dilution factor for production in the MD A and RD' eras. 

Combining Eq. (25) — Eq. (29) gives 


Y 0 i (T) 


A msMpi ( -jr 

jr 
,m B 


T 2 l 



3 

X" ’ 


TrT 1 / 2 


T 


3/2 


M 



(30) 
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For decays and scatterings in all eras, FI via interactions of dimension 4 is IR dominated at 
T ~ rn B . Decays and scatterings then give contributions with the same scaling behavior, giving 
a total yield today of 


Y 


X 2 M J 


PI 


0 i 


m B 


n 7 ’ 
B 


1, E -f- 


T R m J/ 2 


m 


T, 


3/2 


M 



(31) 


Including numerical factors the decay contribution is larger than the scattering one, which is 
why we have included only decays in our numerical analysis. 

For FI during the MD^a era, Eq. (30) shows that IR domination is extreme, with a power 
behavior of 1/T P with p = 8 or 10, so that FI production is maximized at temperature m B /x R j 
with xfi ~ 4. The high power dependence leads to an abundance from FI on the exponential 
tale of the Boltzmann distribution for B that is enhanced by a large factor E, of order 10 4 , 
compared to FI at T ~ m B . Similarly, in Fig 5 Y X (T) peaks at x > 1. The Tj { behavior for the 
yield during the MDtva era, shown in Eq. (31), leads to the straight lines of slope 7 in Fig. 6. 

For FI during the MD^ era, Eq. (31) shows that the yield is proportional to T R /T^ 2 . so 
that in the MD^ region of Fig. 8 the contours of constant t b have D 3 T R constant. On the other 
hand, in the MD^va region of Fig. 8 the contours are vertical since they depend only on T R and 
not on Tm■ The boundary between the MD^a and MDa eras is given by 

(32) 


The boundary between the MD^ and RD' eras is given by rn B rv_/ t m — dt r . 

For any ( Tm,T r ), knowing which era i dominates FI allows the contours of Fig. 8 to be 
normalized by requiring the observed DM abundance 


Y 0i m x ~ T eq ~ leV, 

so that the B decay rate prediction scales as 


Bi 


Xfrri B 


T eq m 


1, 


m 


B 


T] 


-3/2 


M 


ETP 


Trui 1 / 2 ’ 



Mpi nix 

Including relevant numerical factors, these scaling law predictions can be written as 


ct b 

m. 

ct b 

m. 

ct b 

m. 

ct b 

m. 


3 x 10 6 
10 


300 GeYV 


T, 


m B 

7 


mx 


100 GeV 


R 


10 GeV 


300 GeV 
•m B 


m x 


Tj 


100 GeV 

r \ /10 5 GeV\ 3/2 / 300 GeV \ 3/2 / rn x 
eV / \ Tm / V m B / V100 Gi 


GeV, 
T r 

10 MeV 


M 

30 GeV 


T, 


M 


rn B 
3 TeV 
rn B 


GeV 


m x \ 
10 GeV ) 


(RD) 

(MDm) 

(MDA 

(RD') 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 
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leading to the behaviors of Figs 6, 8 and 9. The value of the coupling A* required for dark matter 
genesis during era i can be estimated from (34). 

For FI during era i to give the observed DM abundance, the peak abundance reached before 
dilution is 


yPeak 
* i 


T 

- 1 eq 

m x 


A 


(39) 


with A evaluated at T ~ ms/xFi fc> r the era. This result for Y Peak breaks down for 

MData, MD^, RD' when it gives a result larger than the equilibrium abundance, Y eq , so that 
FI is no longer the appropriate description of the DM production mechanism. This happens for 


(MDrf.RD') 

+ R J-eq 


m B > (yeq rnx] 1/5 

%fiTr ~ \ T eq J 


(MDtva) 


(40) 


and corresponds to the shaded regions labelled “No Freeze-in” in Figs. 8 and 9. Substituting 
these bounds into (34) one discovers that the couplings while much larger than \rd, are 

small. 


3.4 Freeze-In from higher dimensional operators 


As the dimension d of the interaction inducing B —> X increases, so the scattering cross section 
for this transition grows more rapidly with temperature. This has the tendency to convert 
FI from an IR process dominated by temperatures near mg to a UV one, dominated at some 
temperature Tuy rag. For example, a d = 5 operator with coefficient 1/M inducing FI during 
era i leads to a total yield today of 


Y r 


0; 


Mpi Tjjv 
M 2 


I; A 


1r 


T l r 7 “ll /2 rji 

R J-UV 1 R 


m 


6 T ’ 
B 1 UV 


T 


-3/2 


M 


Tj 


M 


(41) 


replacing the result (31) for d — 4. However, in the case of the MDjv^ era, the IR dominance for 
d = 4 is so strong that the scattering remains IR dominated for all d < 8 [5]. The FI result is 
obtained from the d = 4 case by replacing A with mp/M. The feeble strength of the interaction 
necessary for dark matter to originate from FI can be understood as arising from a high mass 
scale M > mg of size 


M 2 ~ ^ T R m x 
T eq m% 

The resulting lifetime, Tp = Tb(Tr, mj^ng), has the same form as for d = 4 shown in Fig. 6, 
except the normalization is changed because FI may be dominated by scattering rather than 
decays and as d increases the enhancement factors A decrease. 



4 Supersymmetric Theories with Axino Dark Matter 

The QCD axion is the most elegant solution to the well known CP problem of strong interactions. 
If supersymmetry is relevant to our universe the details of how the axion solution is implemented 
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can play a crucial role in electroweak physics. This happens if the HuHd operator is charged 
under the Peccei-Quinn (PQ) symmetry whose breaking delivers the QCD axiom If this is the 
case the usual supersymmetric /i-terrn is not allowed in the superpotential but must be generated 
after PQ breaking by operators of the form [20] 


W D 


O^HuHp 

Mr 1 


(43) 


where O ^ is a dimension n operator which we assume to be the product of SM gauge singlet 
superfields with the right PQ charges. The effective /i-terrn comes from the vev of O hd 




(OW) 
Mr 1 ' 


(44) 


The fields that break the PQ symmetry have PQ charges qi and vevs V) and V 2 = of 
is the effective PQ breaking scale. The axion superfield A 


A= S -^ + V20a + 0 2 F (45) 

v2 

transforms under a PQ transformation with real parameter a as A —> A + iaV. The axion 
supersymmetric effective field theory just below V has self interactions and interactions with 
the Higgs supermultiplets given by [21] 


C(A,H UtD ) = I d A d^V 2 e q ' {A+A])/v + J d 2 9e qoA ' v H u H D + h.c. (46) 

where qo is the PQ charge of the operator 0^ n \ Expanding this Lagrangian gives 

Jd 4 e (aa* + ^t + h.c. + ..}j+ Jd 2 eliHuHo (i + ^ + f^r + ■••)+ h.c. (47) 

where b = 'Yh i q^V 2 /V 2 . Additional terms in the Lagrangian give interactions of the axion 
supermultiplet with supersymmetric gauge field strengths W a 

C(A,W a ) = f d^Y.-Caj^ AW aa W a a . (48) 

a 

To compute axino FI we work in a basis where the coefficients c a get contributions only from 
integrating out matter, charged under both the SM gauge group and PQ symmetry, that has a 
mass larger than the scale V. 
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4.1 Axino production in standard cosmology 

In the limit of of unbroken SUSY, the supermultiplet in Eq. (45) contains degenerate particles, 
that are massless up to a small non-pertubative QCD contribution. The axino a acquires a highly 
model dependent mass after SUSY breaking. Given its very weak interactions with visible fields, 
the axino could easily have dramatic consequences for cosmology regardless of where it sits in 
the spectrum of superpartners [22] . We take the axino to be the stable Lightest Supersymmetric 
Particle (LSP), compute its cosmological abundance, and discuss constraints on the theory that 
allow it to be dark matter. In this sub-section we assume a standard cosmological evolution, i.e. 
a non-adiabatic matter dominated era following inflation starting at a temperature T max and 
ending at a temperature Tr, followed by radiation domination down to the eV era. 

Axino production in the early universe can take place through different mechanisms. Out-of- 
equilibrium LOSP decays naturally provide cold axino dark matter [23]. Scattering and decays 
of particles in the primordial plasma give additional contributions to the axino relic density. 
This can be due to processes mediated by the super-gauge axino-gluino-gluon interaction [24] 
or the axino-higgsino-Higgs Yukawa coupling in DFS models [25]. 

Assuming T R > /x, Eqs. (47) and (19) allow us to calculate the axino FI yield coming 
from higgsino decay h —s- ha. In particular the relevant dimension 4 operator in (47) gives 
A = (qo/ 2)(/i/U). With n mg, ms, and negligible mixing with the electroweak gauginos, the 
decay widths of the charged and neutral Dirac components are equal and given by 

r(/A W ±a) = r (h° -+ Z~a , ha) * (49) 

The decay widths do not depend on tan f3 under the assumption of unmixed higgsinos. These 
widths define the higgsino lifetime by r = r -1 . Notice, furthermore, that if the B particle is the 
higgsino g R = 8 in Eq. (19), giving 


>a(FI) 


«9x 10 ~ 4 q 2 0 


/106.75 
V g* 


3/2 


h 


300 GeV 


^10 10 GeV^ 2 


(50) 


An additional contribution to the yield of a arises if some of the coefficient c a , say C 3 , are non¬ 
vanishing. This contribution is UV dominated and arises from processes of the kind GG —> Ga 
where G is either a gluon or a gluino and has been calculated in [26] 


Y a (UV) « 1.7 x 10" 3 c 2 


2 x 106.75 

g* 


3/2 


- R 


10 5 GeV 


10 10 GeV 
V 


(51) 


This growth with T r stops once Tr exceeds the mass of the heaviest fermion carrying both SM 
and PQ charges [27], so that this UV contribution is highly model dependent. In contrast, the 
FI yield (50) is generic to all DFS theories. The axino abundance is given by the FI mechanism 
only if Uj(FI) > Ya(UV) and W(FI) < Y a eq « 2 x 1CT 3 , the thermal abundance . 6 * * 

6 An additional UV dominated contribution to the axino abundance comes from dimension 5 operator of the 

form r/V 2 aaHuHr). These contribution are suppressed by a factor g 2 /V 2 so that they are irrelevant for weak 

scale supersymmetry for the values of Tr we are going to consider. 
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Defining the axion decay constant, /, by the interaction with the gluon held strength G in 
the low energy effective theory, / 32tt 2 ) (a / f)GG, gives 

1 _ c 3 + cj R + 3 qp 

f V2V ’ ( } 


where c^ R is the contribution to the anomaly from non SM fields whose masses are below the 
scale / and the last term in the sum accounts for the anomalous contribution coming from 
the SM quarks. For the choice of parameters made in Eqs. (50) and (51), the axino yield 
is very high. The observed dark matter density will be exceeded unless the axino mass is 
far below the supersymmetry breaking scale and/or / is very large. This is especially true if 
T r > 2.10 5 GeV/cg, when the yield is dominated by UV scattering. Large / and/or small 
are certainly possible, and are the conventional approach to this axino over-abundance problem, 
but they place severe restrictions on the theory. 

The simplest theories giving the superpotential as in Eq. (43) have O O with dimension 
n = 2, giving / 2 < fiMpi for M* < M PP so that the yields of Eqs. (50) and (51) are indeed 
very large. For theories with larger n, a large value of / yields too much axion dark matter 
unless the initial axion misalignment angle is adequately tuned to small values. The tree-level 
axino mass is bounded from above by < m 2 /f [28], with m ~ /i the scale of SUSY breaking. 
This limit applies even taking into account the soft SUSY breaking terms [22]. This mass would 
be small enough to avoid any overabundance. However, this tree-level prediction can easily be 
spoiled. It is quite generic for the axino to pick up a mass of order the gravitino mass via higher 
dimensional operators coupling the axion supermultiplet to SUSY breaking fields [29]. If SUSY 
is broken at a high scale, as in gravity mediation, the axino will generically pick up a weak scale 
mass. Even for lower mediation scales there can be loop-mediated SUSY breaking contributions 
to the axino mass (first pointed out in Ref. [30] within the context of supergravity theories) 
depending on the UV completion of Eq. (47). 

In this section we pursue an alternative solution to this axino-overproduction problem: we 
take Tr below the higgsino mass so that FI occurs during the MD era. The FI yield (50) is then 
suppressed by a combination of less available time for production and subsequent dilution from 
entropy production, as in (31). From the above discussion we assume the axino mass is not far 
below the weak scale 


HT 3 < — < 1. 


(53) 


A 


The UV contribution (51) is similarly suppressed and, as discussed in the next sub-section, 
we will be interested in situations where it is sub-dominant. We find that axino dark matter 
then predicts a restricted range of / and that in much of the parameter space the reactions 
h —> ha , Za and —>• W ± a lead to displaced signals at colliders. 


4.2 Axino Freeze-In during an early matter dominated era 

Let us review the assumptions which will go into this section. We assume the strong CP problem 
to be solved by a supersymmetric DFS axion. The axino mass is required to be in the range 
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of Eq. (53) and we assume the axino to be the DM. We furthermore assume for simplicity the 
LOSP to be a higgsino with negligible mixing with the electroweak gauginos. 

We begin by requiring there are no fields charged under both the PQ and SM gauge sym¬ 
metries that are heavier than the TeV scale, so that c a , = 0 in Eq. (48) and there is no 
UV scattering yield from temperatures above the TeV scale. In this case, after the rescaling of 
Eq. (52), the higgsino decay widths depend only on /, 

r(h ± W ± a) = r(h° Za, ha) « (54) 


The general analysis of FI in a cosmology with an intermediate MD era performed in Sec. (3) 
can then be directly translated to this model using Eq. (54), for both the numerical and analytic 
analyses. For any point of the general parameter space depicted in Fig. (2) we are thus able to 
determine the value of V, and hence /, that gives the observed dark matter abundance. 

The numerical results for axino FI are shown in Figs. (10) and (11) for FI during inflationary 
reheating and during a generic MD cosmology respectively. The qualitative behavior of the 
figures can be understood using the analytical approximations in Sec. (3.3), in particular (31), 
and taking Qb — 8 for the higgsino: 
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(MDjy^) 

(MD a ) (55) 
(RD ; ) 


In the region to the right of the black dashed line in Figs. (10) and (11), there is an over¬ 
abundance of QCD axions produced by the misalignment mechanism. This limit applies to the 
region of parameter space with a post-inflationary scenario [31, 32] where the PQ phase tran¬ 
sition occurs after inflation. If Tr is above 160 MeV the entropy generated at the end of the 
MD era does not dilute the axion density coming from coherent oscillations of the axion held 
but only the (theoretically uncertain) contribution coming from the decay of topological defects 
[33] . For values of Tr below 160 MeV the misalignment contribution is strongly diluted and 
negligible. Notice that the red exclusion does not apply for the pre-inhationary axion scenario 
as one can tune the initial misalignment angle to arbitrarily small values. 

The yellow shaded region in the figures corresponds to the bound coming from axion emission 
from white dwarfs [34], The bound applies if the axion-electron coupling 

Yj d ^ a (56) 

is present, which is typically the case for DFS axions, and is shown for C e = 1/3. The axion 
energy density bound from misalignment and the white dwarf bound depend on / and are shown 
in Figs. (10) and (11). 
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jj = 300 GeV 



li = 3 TeV 



Figure 10: / as a function of Tr that yields axino dark matter for ji = 300 GeV (left) and 
/i = 3 TeV (right). In each case we show results for three different axino masses. Above 
the black dashed line the axion abundance is above the observed value in the post-inflationary 
case, while the dot-dashed lines refer to the pre-inffationary scenario for various choices of the 
initial misalignment angle. The yellow region is excluded by axion emission from white dwarfs. 
Displaced collider signals occur in the shaded horizontal bands as described in Table 1. 


The green shaded regions correspond to higgsino lifetimes that give rise to displaced signals 
at colliders (with the color code as defined in Tab. 1). The higgsino decays to missing energy 
and, depending on whether it is charged or neutral, to a W boson or democratically to a Higgs 
or Z boson. An early study of these signals can be found in Ref. [35], a more recent analysis for 
LHC collisions was performed in Ref. [36]. Notice that this signal is distinguishable from the 
decay of the neutralino NLSP in gauge mediated models where final states containing a photon 
are possible while they are strongly suppressed in the axino case. 

At a hadron collider, in the most favorable situation, it will be possible to measure both ji and 
nia from kinematical distributions and the higgsino decay width from the displaced vertex. This 
is enough to extract Tr in the case of FI during post-inflationary cosmology or a corresponding 
curve in the ( Tr , D) plane in the case of a generic MD era. In the axino case the decay width is 
related to the axion decay constant / and this correspondence is 1 to 1 for an unmixed higgsino 
as in Eq. (54). Even though the axion is not dark matter in most of the displaced vertex 
region, low energy experiments looking for axion-mediated long-range forces [37] can measure / 
independently and provide a crucial verification of the framework discussed in this section. 

How are these results affected by going to non-minimal DFS models with matter heavier 
than the TeV scale inducing c a , c JR ^ 0 in (48) and (51)? For the case of the MD W i cosmology, 
as occurs after inflation as shown in Fig. (3), there is no change and Fig. (10) applies unaltered. 
In this case, even though the operator of (48) has dimension 5, it induces an axino yield that is 
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Figure 11: Contours of / that give axino DM via FI in the (Tr, D) plane for different values of /i 
and rrig. To the right of the black dashed line the axion abundance is above the observed value 
in the post-inflationary case, while the dot-dashed lines refer to the pre-inflationary scenario for 
various choices of the initial misalignment angle. The yellow region is excluded by white dwarfs 
constraints. Displaced collider signals occur in the green and blue shaded bands, as described 
in Table 1. 
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IR rather than UV dominated, and this IR contribution is always less than that from higgsino 
decay. Furthermore, additional FI yields from heavy colored states are suppressed by the strong 
IR dominance of MD^a- However, gauge scattering axino production is UV dominated during 
MD 4 and RD' eras, as well as during any RD era that preceded them, and further FI yields 
from these eras might be significant. In this case it is model-dependent whether the axino yield 
is dominated by FI from higgsino decays and Fig. (11) applies only if the FI contribution from 
higgsinos dominates. 

Figs. (10) and (11) demonstrate that, in a wide range of DFS models, the dark matter can be 
axinos with a mass in the GeV to TeV range, with the dominant production mechanism being 
FI from higgsino decay during an early MD era. Over a large fraction of the parameter space 
not excluded by axion overproduction and white dwarf constraints, there are displace signals at 
colliders arising from precisely the decays that induced the cosmological abundance, h —» ha, Za 
and h± —> W ± a . The relevant range of / is low: from the experimental lower bound up to 
order 10 12 GeV, which can be independently explored by a variety of techniques [37]. 
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A Solution of FI Boltzmann Eq. for Any History 


We solve the Freeze-In Boltzmann equation (17) by employing techniques analogous to the ones 
used in Refs. [4, 5]. In cases where entropy is not conserved in the evolution, Y x is not the best 
variable to use. We define the X comoving number density, X = nxct 3 , and we use the scale 
factor as time variable dt = H~ 1 dlna. 

The Freeze-In Boltzmann equation in terms of the new variable reads 

dx _ r w n ( T ( a )) a3 K \[m B /T{a)\ 
din a B H(a) K^m B lT[a )] 

The asymptotic comoving density is found by imposing the initial condition nx(0) ~ 0 


x 00 = r 


B 


d In a 


, n g(T(a'))a' 3 K 1 [m B /T(a')\ 
H(a') K 2 [m B /T{a ')] 


(58) 


The Hubble parameter H(a') is given by the Friedmann equation 


( 8vr\ 1/2 ypR(a') + pNR(a') 

\ 3 / Mpi 


(59) 
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where pxn is the non-radiation component of the energy density and the temperature is related 
to the radiation energy density pn by Eq. (3). 

For late times the comoving density X approaches a constant value and this value can 
be used to compute the DM density today. At Tr the universe enters the usual RD era and we 
can trust entropy conservation again. Therefore, we identify a late temperature less than Tr 
and evaluate px and the entropy sa 3 . This ratio in the subsequent evolution of the universe 
must be constant. We define 


PxiT < Tr) _ nx(T < Tr) _ X(T < Tr) _ X^ 

s{T < Tr) mX s(T < Tr) mX S(T < Tr) ^ S(T < Tr) ' 


(60) 


The energy density today reads 


Px(to) = (x s 0 , So = 2891 cm 


(61) 


The contribution to the dark matter energy density reads 

fl x h 2 = , Pcr/h 2 = 1.054 x 1CT 5 GeVcm -3 . (62) 

Pci/h 2 

This solution is completely general and valid for Freeze-In during an arbitrary background 
of radiation and non-radiation components. 


B Irrelevance of LOSP Freeze-Out in MD Era 


So far we have not discussed the LOSP Freeze-Out (FO) contribution to the DM density. If we 
assume the LOSP particle B to be a generic particle with a weak scale mass and couplings, and 
in the absence of dilution, this contribution is in the right ballpark to account for the observed 
DM abundance. The goal of this Section is to show that, in the displaced signal regions for FI 
identified above, FO is always suppressed and can be safely neglected. 

The LOSP FO calculation can be performed following the strategy of Refs. [4, 5]. The 
Boltzmann equation describing the evolution of the number of B particles reads 

d ^~ + 3 Hn B = -(<7u r ei) B B—>aii [n% - n e £ 2 ] . (63) 

We parameterize the annihilation cross section as follows 


(<Tn r el)sB-s.all — 
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2 

B 


(64) 


As shown in App. A, when there is dilution the variable Yr is not convenient, since entropy 
is not conserved. Analogously to the previous case we define B = n B a 3 and the Boltzmann 
equation (63) in terms of the new variable B and with time evolution accounted for by the scale 
factor results in 


dB 
din a 


Wei) 

H(a)a 3 


[B 2 {a)-n^{T{a)) 2 a & ] , 


(65) 
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Figure 12: Suppression of LOSP FO contribution to the X relic density as a function of T R for 
inflationary reheating. We show two different choices of ( ms,Tnx ) masses. Displaced collider 
signals occur in the shaded horizontal bands as described in Table 1. 


with the Hubble parameter as in Eq. (59). For large values of the scale factor, corresponding 
to temperature much greater than T R , the variable B approaches an asymptotic value The 
present energy density of A" particles is computed analogously to the Freeze-In case 

(66) 

Notice the difference with the Freeze-In case in Eq. (60) where B^ is replaced by TV,. 

In what follows we consider different cosmologies and different values of m R and mx, and 
for each case we set a R in such a way that in the absence of dilution Eq. (66) reproduces the 
observed DM density. We find the following values 


f 0.07 ( m B ,mx ) = (300 GeV, 100 GeV) 

\ 0.7 (iri8,%) = (3 TeV, 1 TeV) 


(67) 


Then we explore the impact of dilution on these reference points. 

We show our results for the case of inflationary reheating in Fig. 12. The red curve gives the 
dilution factor as a function of Tr over the range of T R where FI can occur, with the two panels 
corresponding to the two cases of Eq. (67). The shaded regions denote the associated collider 
signal for the chosen masses and for the value of T R on the red curve, as found in Fig. 6. 

Next we consider the more general case of both MD a and MD^ a, with results in Fig. 13 for 
the same choices of B and X masses. The contours of both panels show that the FO contribution 
is always reduced as either Tr is reduced or D is increased. The shaded regions of both panels 
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Figure 13: Suppression of LOSP FO contribution to the A" relic density in the (Tr, D) plane 
for two different choices of for lighter X the suppression will be greater. Displaced 

collider signals occur in the green and blue shaded bands, as described in Table 1. 


refer to the FI contribution, and are taken from the appropriate panels of Figs. 8 and 9. In 
the green regions, where FI gives displaced vertices, the FO contribution is suppressed by 1 to 
6 orders of magnitude compared to conventional FO during the RD era, and hence is typically 
expected to be negligible. 

In the right panel of Fig. 13, the three dashed red lines separate the regions where FO occurs, 
from right to left, in the (RD, MDjv. 4 , MD^, RD') eras. These three boundaries, from right to 
left, correspond to FO temperatures Tp of Tr,Tna and Tm■ The FO abundance in region i is 


Y i (, n 

VO ‘ Mpi (av) Tp ( T|’ 

which can be compared to the FI result of Eq. (31) 
contours in both panels of Figs. 12 and 13. 


Tr T r \ 

Ttfrr Tm) 

and gives the asymptotic slopes of the 



References 

[1] S. Hannestad, Phys. Rev. D 70 , 043506 (2004) [astro-ph/0403291]; F. De Bernardis, Nucl. 
Phys. Proc. Suppl. 194 , 39 (2009). 

[2] M. Kawasaki, K. Kohri and N. Sugiyama, Phys. Rev. D 62 , 023506 (2000) [astro- 
ph/0002127]; K. Ichikawa, M. Kawasaki and F. Takahashi, Phys. Rev. D 72 , 043522 (2005) 
[astro-ph /0505395]. 


29 
































[3] J. McDonald, Phys. Rev. D 43, 1063 (1991). 

[4] D. J. H. Chung, E. W. Kolb and A. Riotto, Phys. Rev. D 60, 063504 (1999) [hep-ph/9809453]. 

[5] G. F. Giudice, E. W. Kolb and A. Riotto, Phys. Rev. D 64, 023508 (2001) [hep-ph/0005123]. 

[6] M. Kamionkowski and M. S. Turner, Phys. Rev. D 42, 3310 (1990). 

[7] L. J. Hall, K. Jedamzik, J. March-Russell and S. M. West, JHEP 1003, 080 (2010) 
[arXiv:0911.1120 [hep-ph]]. 

[8] T. Moroi and L. Randall, Nucl. Phys. B 570, 455 (2000) [hep-ph/9906527]; R. Allahverdi and 
M. Drees, Phys. Rev. Lett. 89, 091302 (2002) [hep-ph/0203118]; R. Allahverdi and M. Drees, 
Phys. Rev. D 66, 063513 (2002) [hep-ph/0205246]. B. S. Acharya, G. Kane, S. Watson and 
P. Kumar, Phys. Rev. D 80, 083529 (2009) [arXiv:0908.2430], 

[9] A. Monteux and C. S. Shin, arXiv: 1505.03149 [hep-ph], 

[10] C. Cheung, G. Elor, L. J. Hall and P. Kumar, JHEP 1103, 042 (2011) [arXiv: 1010.0022]. 

[11] S. Chatrchyan et al. [CMS Collaboration], JHEP 1302, 085 (2013) [arXiv:1211.2472 [hep- 
ex]]. 

[12] G. Aad et al. [ATLAS Collaboration], arXiv: 1504.05162 [hep-ex], 

[13] V. Khachatryan et al. [CMS Collaboration], Phys. Rev. D 91, no. 1, 012007 (2015) 
[arXiv:1411.6530 [hep-ex]]. 

[14] G. Aad et al. [ATLAS Collaboration], arXiv: 1504.03634 [hep-ex], 

[15] G. Aad et al. [ATLAS Collaboration], JHEP 1411, 088 (2014) [arXiv:1409.0746 [hep-ex]]. 

[16] V. Khachatryan et al. [CMS Collaboration], Phys. Rev. D 91, no. 5, 052012 (2015) 
[arXiv:1411.6977 [hep-ex]]. 

[17] G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 88, no. 11, 112003 (2013) 
[arXiv:1310.6584 [hep-ex]]. 

[18] V. Khachatryan et al. [CMS Collaboration], Eur. Phys. J. C 75, no. 4, 151 (2015) 
[arXiv: 1501.05603 [hep-ex]]. 

[19] S. Dimopoulos and L. J. Hall, Phys. Lett. B 196, 135 (1987). 

[20] J. E. Kim and H. P. Nilles, Phys. Lett. B 138, 150 (1984). 

[21] B. Bellazzini, C. Csaki, J. Hubisz, J. Shao and P. Tanedo, JHEP 1109, 035 (2011) 
[arXiv:1106.2162 [hep-ph]]. 

[22] K. Rajagopal, M. S. Turner and F. Wilczek, Nucl. Phys. B 358, 447 (1991). 


30 



[23] L. Covi, J. E. Kim and L. Roszkowski, Phys. Rev. Lett. 82, 4180 (1999) [hep-ph/9905212], 

[24] L. Covi, H. B. Kim, J. E. Kim and L. Roszkowski, JHEP 0105, 033 (2001) [hep- 
ph/0101009], 

[25] E. J. Chun, Phys. Rev. D 84, 043509 (2011) [arXiv: 1104.2219 [hep-ph]]. 

[26] A. Strumia, JHEP 1006, 036 (2010) [arXiv: 1003.5847 [hep-ph]]. 

[27] K. J. Bae, K. Choi and S. H. Im, JHEP 1108, 065 (2011) [arXiv: 1106.2452 [hep-ph]]. 

[28] K. Tamvakis and D. Wyler, Phys. Lett. B 112, 451 (1982). 

[29] C. Cheung, G. Elor and L. J. Hall, Phys. Rev. D 85, 015008 (2012) [arXiv: 1104.0692 
[hep-ph]]. 

[30] P. Moxhay and K. Yamamoto, Phys. Lett. B 151, 363 (1985). 

[31] M. P. Hertzberg, M. Tegmark and F. Wilczek, Phys. Rev. D 78, 083507 (2008) 
[arXiv:0807.1726 [astro-ph]]. 

[32] F. D’Eramo, L. J. Hall and D. Pappadopulo, JHEP 1411, 108 (2014) [arXiv:1409.5123 
[hep-ph]]. 

[33] R. L. Davis, Phys. Lett. B 180, 225 (1986); C. Hagmann, S. Chang and P. Sikivie, Phys. 
Rev. D 63, 125018 (2001) [hep-ph/0012361]; T. Hiramatsu, M. Kawasaki, T. Sekiguchi, 
M. Yamagnchi and J. Yokoyama, Phys. Rev. D 83, 123531 (2011) [arXiv:1012.5502 [hep- 
ph]]; T. Hiramatsu, M. Kawasaki, K. Saikawa and T. Sekiguchi, Phys. Rev. D 85, 105020 
(2012) [Erratum-ibid. D 86, 089902 (2012)] [arXiv:1202.5851 [hep-ph]]. 

[34] K. A. Olive et al. [Particle Data Group Collaboration], Chin. Phys. C 38, 090001 (2014). 

[35] S. P. Martin, Phys. Rev. D 62, 095008 (2000) [hep-ph/0005116], 

[36] G. Barenboim, E. J. Chun, S. Jung and W. I. Park, Phys. Rev. D 90, no. 3, 035020 (2014) 
[arXiv: 1407.1218 [hep-ph]]. 

[37] A. Arvanitaki and A. A. Geraci, Phys. Rev. Lett. 113, no. 16, 161801 (2014) 
[arXiv:1403.1290 [hep-ph]]. 


31 



